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A linear increase with the volume of the topological susceptibility can signal spontaneous breaking of parity P 
and time inversion T, due to a nonzero vacuum expectation value of the topological charge Q. Such a breaking 
would produce a nonzero electric dipole moment of the neutron, d n . An upper limit to d„ is derived from numerical 
simulations at increasing volumes. 



1. INTRODUCTION 

The QCD lagrangian 

£qcd = -TriG^Gfiv} + 

^2i} f (ip-m f )ipf (1) 
/ 

is invariant under parity P and time inversion 
transformations T. Adding to £qcd a Chern- 
Simons term 

£ cs = 0Q(x) , (2) 
with Q(x) the topological charge density 

g 2 

Q( x ) — 2 e ^vp<j^{G tiy G pa } (3) 

does not modify the equations of motion since 
Q(x) is a divergence, Q(x) — d t _ l K tJ/ (x), 

K„ = (d P Al - l -gf abc A b p A^ . (4) 

The lagrangian Ccs has a dynamical effect, be- 
ing sensitive to the boundary conditions, and its 
presence in C = £qcd + £cs results in an explicit 
breaking of P and T. Indeed, since Q ~ E ■ H it 
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is odd under both of them. This breaking can 
produce a nonzero value of d n , the electric dipole 
moment of the neutron [J, d n « 8 ■ lO" 16 ^ ■ cm. 
The experimental upper limit is cZ n | exp < 

6 ■ 10 _26 e ■ cm which implies (9~1CP 10 . This is 
a small number and its smallness is not protected 
by any known symmetry. 

An alternative possibility for having d n ^ 
could be a spontaneous breaking of P by a nonzero 
value (0|Q|0). This possibility however would be 
excluded by a theorem by Vafa and Witten [3]. 
Nevertheless this theorem is based on assump- 
tions whose validity has been questioned. We 
have then addressed the issue of a spontaneous 
breaking of P and T by numerical simulations of 
QCD on the lattice Hj. 

2. THE VAFA WITTEN THEOREM 

The statement is that in a theory of gauge fields 
and fermions like QCD which is Lorentz invariant 
and with vector interactions, the v.e.v. of any op- 
erator O odd under parity is zero. The argument 
goes schematically as follows: let us consider the 
generalized partition function 




where ip represents the generic field variables, and 
its continuation to the euclidean 

Z Euclid [d]=e- VE ^ e \ (6) 

where E{9) is the free energy as a function of the 
9 parameter and V is the spacetime volume. 

After the analytic continuation the measure is 
positive since the action is real and the fermion 
determinant is positive. The operator O instead 
acquires an extra factor i since, to be odd under 
P and T, it must have an odd number of invariant 
tensors €^ vp(J , at least if it depends only on the 
gauge fields like Q(x). 

As a consequence the following equality holds 



-VE{6) _ 



Vip e 



-S(tp)-i$ J Od 4 x 



(7) 



If O is hermitean then the term exp (— iO J Od 4 x) 
becomes a phase factor and since the rest of the 
integrand is positive, the inequality 



-VE(8) < -VE{0) 



(8) 



derives. Consequently we have E(9) > E(0), i.e. 
9 = is a minimum. Since for small 9, E(9) ~ 
1 - 9 (0|Q|0), it implies (0|Q|0) = 0. 

There are a number of assumptions behind this 
result. 

• Lorentz invariance, which is certainly true 
for QCD at zero temperature, is lost at fi- 
nite temperature. Then, the breaking of 
parity, (0|Q|0) ^ becomes plausible |5lfi| 
with possible observable consequences in 
heavy ion collisions. 

• The argument can be nonvalid if O is con- 
structed not only by use of gauge fields but 
also with fermionic fields 0. 

• Analyticity in 9 of Z[9] is assumed. This 
hypothesis has been debated in the litera- 
ture |8l9ll()j and it is in fact the one which 
we have tried to test on the lattice for the 
topological charge operator. 

Suppose indeed that E(9) has a minimum at 
9 = 0, but a discontinuous derivative (a cusp) at 
this point so that 

dE(9) 



dd 



9=o± = ±«, 



(9) 



(0|Q|0}| e=0± =±Va. 



(10) 



Then a spontaneous breaking of P and T would 
result. 

The topological susceptibility 



X = / d A x <0|Q(x)Q(0)|0> = 



(Q|Q 2 |Q) 

V 



(11) 



would then contain the usual contribution coming 
from vacuum fluctuations (which in the present 
paper shall be called x), and another contribu- 
tion derived from the nonzero vacuum expecta- 
tion (0|Q|0), 



x = a V + \ 



(12) 



and x would be infrared divergent. 

The electric dipole moment of the neutron d n 
can be estimated as in pQ. It is proportional 
to the parameter a times a typical volume m~ A 
times a factor (m^jm, n ) which makes it vanish 
in the chiral limit, 



(13) 



An extra factor 1 / m n is needed to assign the cor- 
rect length units. From lattice we will get an 
upper limit for a and, from Eq. an upper 

limit for d n . 

3. LATTICE ANALYSIS 

Lattice is an UV regulator of the theory. The 
topological charge Q can be determined either 
directly for the discretized version Ql of the 
charge operator (the so-called field-theoretical 
method) or from the counting of zero modes in 
the fermionic degrees of freedom by using the 
anomaly equation 



d^h tl = 2N f Q{x). 



(14) 



The latter method requires a special discretiza- 
tion of the fermion fields ^2] which at the present 
state of the art is rather demanding in computer 
time, especially when one has to deal with large 
volumes as in our case. 



In the quenched theory the general relations 
Ql = ZQ , X L = ^L = Z 2 a 4 x + M (15) 

hold. Z is a multiplicative constant C2| that can 
be determined by the expression 

Z=-(Ql)\q=u (16) 
n 

where the v.e.v. is calculated in the topological 
charge sector Q = n where n is a nonzero inte- 
ger. The procedure of selection of the set of con- 
figurations belonging to this sector is performed 
by cooling. On the other hand M is an additive 
renormalization due to singularities that arise af- 
ter the limit x — > in the integrand of Eq. . 
M is an UV effect that can be singled out by cal- 
culating \L in the zero charge topological sector, 

M = X l\q=v EH] 

The Ginsparg- Wilson formalism allows to pre- 
serve exact chiral symmetry in the theory. In 
this case Z = I and M = Q3j|. In the field- 
theoretical method these renormalization con- 
stants must be calculated with the help of cool- 
ing-heating methods [14115) . 

We have adopted the field-theoretical method. 
The simulations have been done on volumes 16 4 , 
32 4 and 48 4 at /3 = 6/g£ are = 6.0, i.e. at lattice 
spacing a w 0.1 Fermi. The statistics was 12 • 10 4 , 
6 • 10 4 and 5 • 10 4 independent measurements for 
the three volumes respectively. 

We have first determined the dependence of Z 
and M on the lattice size L = V 1 ' 4 . In Fi gures 1 
and 2 we show the corresponding results. Only 
the additive renormalization M, which is an ex- 
tensive quantity, displays a clear size dependence. 

Then, knowing Z, M and the lattice spacing a 
we have extracted \ as a function of the volume 
V = L 4 . The result is shown in Figure 3. 

No linear dependence on L 4 is found within 
errors and this allows to put the upper bound 
(within 1-cr) 

a<( — - — 1 , (17) 
~ \4 Fermi/ v ' 

or, by use of Eq. l(T^|) . 

d n < 3.5 10~ 21 e-cm. (18) 
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Figure 1. Z versus L for the lattice definition of 
the topological charge operator Q l adopted in the 
present paper at f3 = 6.0. The line is the result of 
the fit (shown in the legend) and the grey band 
is the 1— a error. 



As a byproduct, we found that the physical topo- 
logical susceptibility at a ~ 0.1 Fermi x ( or X 
if a vanishes) is (l73.4(±0.5)(±1.2)(i^) MeV) 4 
where the errors are the statistical error from our 
data, the one derived from the value used for 
and an estimate of the systematic error respec- 
tively This result agrees with previous determi- 
nations C3 an d with the Witten-Veneziano for- 
mula CHI needed to explain the large i]' mass. 

4. CONCLUDING REMARKS 

Our upper limit on d n is 5 orders of magnitude 
bigger than the experimental one. It can, how- 
ever, be improved by increasing the lattice size 
and by decreasing the statistical errors. 

We are repeating the computation at nonzero 
temperature to check the ideas of reference 
The field-theoretical method to calculate the 
topological susceptibility on the lattice is fast and 
efficient and is presently the only one which can 
afford large lattices. 

The calculation was done on the APEmille fa- 
cility of INFN in Pisa. 

This work is partially financed by MIUR, Pro- 
gram "Frontier Problems in the Theory of Fun- 
damental Interactions" . 
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Figure 2. M versus L for the lattice definition of 
the topological susceptibility xl adopted in the 
present paper at (3 = 6.0. The line is the result of 
the fit (shown in the legend) and the grey band 
is the 1— a error. 
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